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Abstract 

Translocation of a polymer out of curved surface or membrane is studied 
via mean first passage time approach. Membrane curvature gives rise to a con- 
straint on polymer conformation, which effectively drives the polymer to the 
outside of membrane where the available volume of polymer conformational 
fluctuation is larger. Considering a polymer release out of spherical vesicle, 
polymer translocation time r is changed to the scaling behavior r ~ L? for 
R < Rg, from r ~ L? for R S> Rg, where L is the polymer contour length 
and R, Rg are vesicle radius and polymer radius of gyration respectively. Also 
the polymer capture into a spherical budd is studied and possible apparatus 
for easy capture is suggested. 
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I. INTRODUCTION 



Membrane and polymer flexibility is the major intrinsic physical property of biological 
organisms. Typical value of biomembrane bending rigidity k is order of a few ksT which 
determines de Gennes-Taupin persistence length £ ~ a exp(27r K/ksT), where a is the molec- 
ular length scale. Now that £ is usually much larger than molecular scales, conformational 
fluctuation of membrane over length scales smaller than £ can be neglected. As far as the 
membrane effect on proteins concerned, thermal fluctuation of membrane can be neglected 
if the dimension of protein is less than £. 

Subcellular compartments usually have vesicular shapes and they exhibit finite curvature, 
which introduces a new length scale R, the local radius of curvature of membrane. The 
effective interaction between protein and membrane is significantly modified due to the 
membrane curvature. For instance, the adsorption-desorption transition of a polymer on 
curved surface is known to occur at a lower transition temperature It is because the 
entropy of a polymer increases near a convex surface, compared with that near a planar 
surface. In this case, the criteria that determines the importance of membrane curvature is 
given by R < Rg, where Rg is the radius of gyration of the polymer. 

There are three length scales involved(i?G ; R, and £) in membrane-polymer system if 
we consider ideal flexible polymer. In this paper, only the large £ limit (£ ^> R,Rq) will 
be explored to examine the membrane curvature effect on polymer translocation dynamics. 
Then the conformational fluctuation of membrane can be safely neglected and the problem 
is reduced to that with flexible polymer and rigid membrane of nonvanishing curvature. 

We take into account only the steric interaction between polymer and membrane, from 
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which we derive entropic free energy barrier of polymer translocation across curved mem- 
brane. For translocation dynamics, we use our previous model || in terms of Fokker-Planck 
equation, where mean first passage time is obtained as a measure of polymer translocation 
time. As results of our study, we present the effect of chain length and membrane geom- 
etry on translocation time. Two specific examples of different membrane geometries are 
considered to exemplify the curvature effect on translocation dynamics. The first one is 
release of a polymer out of spherical vesicle, and the second one is polymer translocation 
crossing two joined vesicles through a small bottleneck. The main results of this paper are: 
(1) Membrane curvature drives the polymer out of spherical vesicle due to polymer entropy 
effect. (2) Polymer capture into a small budd takes very long time proportional to exp(L), 
which can be reduced to algebraic dependence on L(=polymer contour length) if there is 
segmentwise energetic bias larger than a critical value. 

In section II, equilibrium conformation and entropic free energy of translocating polymer 
are determined as a function of membrane curvature. Dynamics of polymer translocation is 
examined and the mean first passage time is calculated in section III. In section IV, summary 
and conclusion is given. 

II. FREE ENERGY FUNCTION OF POLYMER RELEASE OUT OF A SPHERE 

The free energy of a polymer is determined by the interaction with environments and 
constraints on it. Here we consider the translocating polymer as composed of two indepen- 
dent parts: two end anchored polymers in the outer space of a sphere and in the interior of 
it, respectively(Fig. |TJ). The interaction between the polymer segments is neglected and only 
the impenetrability constraint on the polymer segment is included as a boundary condition 
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of polymer segmental distribution at membrane surface. 

Suppose that there is a polymer located either inside or outside of a sphere. The inside 
polymer is confined in a sphere of radius R, and the outside one is excluded by the same 
sphere. We consider the ideal flexible polymer, which is composed of n segments with Kuhn 
length b either inside or outside a sphere, whose one end position is fixed at r = (r , 9 , 4> ) 
in spherical coordinate. Then the Green's function of this polymer satisfies the following 
Edwards' equation ||, 

G ± (r\r 0] n)=6^\r-r )5(n), (1) 

where r = (r, 9, (f>), and G ± (r|ro; n) denotes the Green's function for the polymer inside and 
outside the sphere respectively. The impenetrability condition due to the sphere can be 
represented by the following boundary conditions, 

G ± (r|r ;n)=0 at r = R. (2) 

Then the partition functions of polymers of n segments located at inner and outer space of 
a sphere are given by 

Z±(r )= / dr r dOsme d^G^fr, r ; n), (3) 
Jn± Jo Jo 

where Q± represents the radial integration range given by (0, R) and (R, oo) for respec- 
tively. Now that physics should be invariant under rotation of the coordinate around the 
origin, is independent of 9q and (f>o, and should be a function of only r . 
Integrating Eq. [1] over 9 and yields the following radial equation, 

Z±(r,r ) = ^6(r-r )5(n). (4) 
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where Z^(r,ro) are the radial Green's functions defined by 
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The solution of this radial equation with the boundary conditions in Eq. is then given by 
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Using these results, we can finally arrive at the partition function of a polymer whose one 
end is fixed at radial position ro as 

Z%(r )= I drZ±(r,r ). (9) 
Jn± 

Let us now consider the partition function of a polymer whose one end is anchored on 
spherical surface. Introducing a sufficiently small anchorage size e, which is used to define 
the anchored end position of the polymer as r = R =F e inside and outside respectively. 
Substituting it into the partition function, we have the following explicit expressions for 
partition functions of end anchored polymers: 
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upto leading order in e. Note that Eqs. O and 11 are the statistical weights of polymers 



anchored on curved surface relative to that in free space. The relative statistical weight only 
due to curvature effect can be obtained as followings: 



which is identical to the result of Hiergeist and Lipowsky [H valid in small curvature limit. 



Using Eqs. [10] and |TTJ, the free energy of the polymer whose one end anchored on surface 
can be obtained as 



F ± {n;R) = -k B T log Z± 

-^TlogEr=iexp(-^A; 2 ) (inside) 

1/21 



(14) 
(15) 



-^Tlog[l + {^} '"j (outside) 
apart from additive constants. Note that these free energy expressions are valid for all 

curvature values. In the limit of R 3> Rq = A^ 1 / 2 6/3, both the inside and the outside free 



energy expressions converge to 



F ± (n;R) 



k B T 



■ log n + constant 



(16) 



which is the conformational free energy of a polymer whose one end anchored on planar 
membrane @. 

For the chain of N segments translocating from the inner side of a sphere to outside, the 



free energy function is given by 

;F( n ) = F-(n; R) + F + (N - n; R) 
= -k B T log 



6R 2 



b 2 



- k B T log X! ex P( 



k=\ 



7i 2 b 2 (N - n) 
6R 2 



k 



(17) 
(18) 



where n is the segment number outside(Fig. |l|). As is depicted in Fig. ^ for different values 
of R, T(n) exhibits nearly symmetric barrier like that of the translocation across planar 
membrane for R ^> Rq, If R becomes comparable or less than Rq, T{n) becomes slanted 
down to the right, which indicate the polymer release is favorable for R < Rq- 
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III. DYNAMICS OF POLYMER TRANSLOCATION 

As shown in [Q, the translocation of a polymer can be thought of as a one-dimensional 
diffusion process of translocation coordinate n, defined by the number of polymer segments 
in the target side, under the effective potential field JF(n). The probability density of n(t) 
given the initial value no is described by Fokker-Planck equation 

d 

— P(n, t\n ) = C FP (n)P(n, t\n ), (19) 

where Cfp(ti) is the Fokker-Planck operator given by 

C FP (n) = __D(n)exp[-#F(7i)]— exp[/3^(n)], (20) 

with D(n) defined as diffusion coefficient. The translocation time of a polymer can be 
defined in terms of mean first passage time r(n; no), time taken for diffusion from no to n, 
which satisfies 

C F p(n )T(n;n ) = -1, (21) 

where C FP (n ) is the backward Fokker-Planck operator defined by 

4pK) = ^exp[/3^(no)]^- J D(n )exp[-^(n )] — . (22) 

Using the following boundary conditions 

^-r(n; n = 1) = 0, (23) 
dn Q 

r(n; no = JV-1) = 0, (24) 

the solution of the above backward equation(Eq. |2l]) can be formally obtained as 

t ee t(N — 1, 1) = ^- / tZne^W / rfn'e-W). (25) 
D Ji Ji 



Here we have set the diffusion coefficient of the whole chain as D(n) = D = ksT ' j '(Nj) ~ 
N^ 1 , where 7 is hydro dynamic friction coefficient for a single segment. Eq. |25| measures 
the diffusion time of the polymer starting from the front segment located in target side and 
ending up only the last segment remaining in the incipient side. The reflecting boundary 
condition at n — 1 means the front segment cannot cross the pore via backward diffusion, 
and hence it is allowed to be located only in the target side. 

As a first example with curved membrane, let us consider polymer release out of spherical 
vesicle as shown in Fig. [j]. The translocation time(Eq. |25|), using the free energy function 
Eq. [18|, is calculated for the case R 3> Rg'- 
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where A = (8vr/15 - 448/225) ~ -0.3156, L = Nb, and a = (n/2)^ 2 ( y R G /R) with 
R G = Nb 2 /3. The a = limit, r ~ (tt 2 /8)L 2 /(2D) ~ L 3 , is just the planar membrane 
translocation time(Fig. ^| A), which is the same as in 0. As R decreases, r decreases be- 
cause the confinement free energy of inside polymer drives the translocation outwards. This 
effect becomes more prominent when R < Rg- In the limit of a — > 00 i.e. R Rg, the 
translocation time is given by 



T 



2D 



(27) 



Note that the leading term scales as r ~ L 2 since a ~ L 1 ' 2 and D ~ L~ x . This reflects 
the fact that confined polymer is squeezed out to the outside, because the confinement costs 
more free energy than being released off. 

The chain length dependence of translocation time is shown in Fig. ^ for different values 
of R. It is remarkable that the translocation time exhibits a crossover from r ~ L 3 to r ~ L 2 
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near R = Rc{Fig. |3| B, C). Also interesting is that this crossover of length scaling behavior 
is the same as that by chemical potential bias studied by us previously [0]. In both cases, 
crossovers are the consequences of membrane asymmetry which adds a linear term to the 
free energy function. 

In fact, as the polymer segment concentration is high enough when R <C Rq, the excluded 
volume effect (EVE) becomes nonnegligible and it is expected to affect the translocation dy- 
namics significantly. A simple scaling argument gives the following confinement free energy 
expression with EVE J7j: 

■RgV /u 



Fin(N) ~ k B T 



N (28) 



R 

where v ~ 3/5 is the swelling exponent of self avoiding polymer in three dimension. Now that 
this free energy is larger than that without EVE, EVE would definitely enhance the outward 
translocation. But the free energy Fi n (N) is proportional to N just like the ideal polymer 
free energy expression for R <C Rg limit. Therefore, the scaling behavior of translocation 
time discussed above will not be changed due to EVE. Only the prefactor is modified to 
reduce the translocation time. 

As a second example, we consider two joined vesicles of different radius Ri and R2 
between which a small bottleneck is opened as shown in Fig. |] (A). In the limiting situation 
that both R 1 and R 2 much larger than R G , the problem is simply reduced to that across 
the planar membrane. The free energy function of polymer translocation for arbitrary Ri 
and R2 is given by 

JF(n) = F + (n;R 2 ) + F + (N -n;^) (29) 
which is depicted in Fig. [| for various values of R\ and R 2 - For Ri ^> R 2 as shown in Fig. 



fl (B), the problem is reduced to polymer capture(delivery) dynamics into a finite size budd 
from planar membrane side. As the capture proceeds, polymer free energy increases, which 
prohibits the capture process and translocation will take longer time than other cases. On 
the other hand, in the opposite limit of R\ <C R2, the problem becomes polymer release 
to planar membrane side out of spherical confinement. This is qualitatively similar with 
the previous example, where the main driving mechanism is confinement free energy of the 
vesicle. 

The translocation time as a function of arbitrary Rx and R2 has the following form: 

= [^1 l\ ( X J gjk exp(-rai(l - y)k 2 ) Eg°=i eM~™ 2 2 yk 2 ) 
' " [d\J Jo y Y.k=ieM-™l(l - x)k 2 )Y%LieM-™lxW) K } 

where a x = (vr/2) 1 / 2 ( J R G / J R 1 ) and a 2 = (7r/2)^ 2 (R G /R 2 ) with R G = Nb 2 /3. In Fig. §, 
translocation time is shown for various values of a± and ol%. For polymers of short length 
such as R G <C Ri,R 2 , the translocating chain does not feel membrane curvature. In this case, 
translocation time is given by r ~ L 3 , the result of planar membrane translocation. If a± ~ 
a 2 for long chain, the effective potential exhibits symmetric barrier and the translocation 
time becomes r ~ L 3 (Fig. ^ B) again. No dramatic driving mechanism can be seen in this 
case, because there is no asymmetry across the bottleneck. 

For ai > «2, translocation time is changed to r ~ L 2 provided that R G > #i(Fig. § 
C). The confinement of R\ radius vesicle squeezes out the polymer in this regime. On the 
other hand, for ot\ C «2, translocation time rapidly increases if R G > i?2- This signifies 
that spontaneous capture by a small budd(i? 2 < Rg) is hard to occur because the free 
energy barrier height increases linearly as the chain length increases(Fig. |5| A), which finally 
results in exponential increase of translocation time as chain length increases(Fig. ^| A). To 
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overcome this difficulty of capture process, segmental energetic bias or its fluctuation || 
can be an apparatus to make the capture be accelerated. Let us here consider the segmental 
chemical potential difference A/x between the two sides, which will add a new contribution 
on free energy function in Eq. as 

AjF(n) = nApL (31) 

which is identical to that introduced in [0. For A/x < 0, the capture process will be 
accelerated and this effect can be dominant over the oppositely directed entropic bias due 
to membrane curvature provided 

|/3A/x| > |/3A/x c | = ^ (32) 

where |/3A/x c | <C 1. If Eq. [32] is fulfilled, the translocation time scales as r ~ L 3 , or even 
t ~ L 2 for |/3A/x| > |/3A/x c | + 1/N which is also very small. In addition, as is shown in || 
for translocation across planar membrane, chemical potential fluctuation can also enhance 
the translocation dramatically. These signify that minute segmental chemical potential bias 
or its fluctuation can be a nice apparatus to make the capture occur easily. 

IV. SUMMARY AND CONCLUSION 

Membrane curvature effect on polymer translocation is explored within our diffusive 
transport model. The geometrical constraint is found to be important to determine the 
entropic free energy of translocating polymer, and asymmetries given by the membrane 
curvature is found to be a possible driving mechanism of polymer translocation. There occurs 
crossover in chain length dependence of translocation time depending upon the membrane 
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curvature. Further, the excluded volume effect is found to be irrelevant to the chain length 
scaling behavior of translocation time, although it affects the polymer conformation and 
translocation dynamics significantly. Finally, entropically-prohibited polymer capture into 
a budd can be accomplished if a minute chemical potential bias is introduced. 
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FIGURES 

FIG. 1. Polymer release out of spherical vesicle of radius R. The translocating polymer can be 
identified with two end anchored polymers composed of n and N — n segments outside and inside 
of vesicle, respectively. 

FIG. 2. Free energy function ^(n), in units of k B T, of polymer release out of sphere as a 
function of translocation coordinate n. (N = 1000, A. R = WRg, B. R = Rg, C. R = 0.5Rq-) 

FIG. 3. Time r for polymer release (translocation) out of sphere, in units of b 2 /(2Do) with 
D = k B T/-f, versus chain length N. (A. R = 3006, B. R = 306, C. R = 156.) Crossover from 
r ~ L 3 to r ~ L 2 occurs near the N corresponding to R = Rg- 

FIG. 4. Polymer transfer between two joining vesicles with different radii Ri and i?2- Translo- 
cation is considered from R\ vesicle to R2 vesicle. (A. R\ is larger than R2, B. R\ » R2) 

FIG. 5. Free energy .F(n), in units of ksT, of a polymer translocating two spheres as a function 
of translocation coordinate n. (N = 1000, A. a\ = 0.5 and 02 = 2, B. a\ = 02 = 1, C. a.\ = 2 and 
a 2 = 0.5.) 

FIG. 6. Translocation time r, in units of 6 2 /(2Z?o) with Do = ksT/j, between two spheres 
versus chain length N. (R 2 = 306, A. R ± = 606, B. R x = 306, C. R x = 156.) 
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Fig. 1, Park, Polymer Release . 



Fig. 2, Park, Polymer Release 




Fig. 3, Park, Polymer Release 





Fig. 4, Park, Polymer Release 



Fig. 5, Park, Polymer Release 




Fig. 6, Park, Polymer Release 




